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(57) Abstract: The invention relates to a system for the 
automatic determination of the mass density of an object 
(100), said system comprising a unit (2) for determining a 
significant dimension (x) of the object (100), a unit (30) for 
determining the intensity (I) of a photonic beam attenuated 
by the passage of the object (100), an acquisition, process- 
ing and analysis unit (200), means (70, 72, 80, 82, 84, 86, 
88) for transporting the object (100), first means (74, 76, 
78) for adjusting the position of the object (100), and sec- 
ond means (90, 92, 94, 96, 98) for adjusting the position of 
the object (100). The invention also relates to a method for 
implementing the inventive system, said method compris- 
ing steps for calibrating components of units 2 and 30, and 
steps for determining the significant dimension of the ob- 
jects (100). Said steps are performed on each object (100) 
from the batch of objects. 

(57) Abrege : Le systeme, pour la determination automa- 
tique de la masse volumique d'un objet (100) comprend : 
une unite (2) de determination d'une dimension significa- 
tive (x) dudit objet (100), une unite (30) de determination 
de l'intensite (I) d'un faisceau photonique, attenuee par la 
traversee dudit objet (100), une unite d'acquisition, de trai- 
tement et d'analyse (200), des moyens de transport 
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(70, 72, 80, 82, 84, 86, 88) de l'objet (100), des premiers moyens d'ajustement (74, 76, 78) de la position de l'objet (100), des deux- 
i ernes moyens d'ajustement (90, 92, 94, 96, 98) de la position de l'objet (100). Le procede de mise en oeuvre du systeme precedent 
comprend des etapes d'etalonnage de composants des unites 2 et 30, et des etapes de determination proprement dites de la dimension 
significative des objets (100), qui sont effectuees sur chaque objet (100) dudit lot d'objets. 
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PROCEDE ET SYSTEME DE DETERMINATION DE LA. MASSE 
VOLUMIQUE ET DES CARACTERISTIQUES DIMENSIONNELLES D«UN 
OB JET, ET APPLICATION AU CONTROLE DES PASTILLES DE 
COMBUSTIBLE NUCLEAIRE EN COURS DE FABRICATION 

5 

DESCRIPTION 

DOMAINE TECHNIQUE 

La presente invention se rapporte au 
domatine des techniques d' analyse non destructives. 

10 L' invention se rapporte plus precisement a 

un procede et un systeme de determination automatique 
de la masse volumique d'objets par attenuation 
photonique en ligne droite et de leurs caracteristiques 
dimonsionnelles . 

15 Elle trouve une application dans le 

controle et le suivi du bon f onctionnement ci' unites de 
fabrication et d'usinage d'objets, comme par exemple 
des pastilles de combustible nucleaire, tel que UOX 
et/ou MOX, et permet en particulier de suivre la 

20 rep product ibilite de fabrication desdits objects, en ce 
qui concerne leur masse volumique. 

Elle peut aussi etre utilisee pour 
determiner des gradients axiaux et radiaux: de masse 
vol-umique faisant ainsi office de tomodensitometre tres 

2 5 precis. 

ETAT DE LA TECHNIQUE ANTE R I EURE 

Des methodes nucleaires non destructives 
actives de determination de la masse volumiqiae ont. deja 
ete elaborees, notamment pour determiner la masse 
30 volumique d 7 echantillons geologiques . Dans Le document 



WO 2005/100952 



PCT/FR2005/000838 



reference : Been, K. , « Nondestructive Soil Bulk 
Density Measurement by X-ray Attenuation », 
Geotechnical Testing Journal, GTJODJ, Vol. 4, No. 4, 
Dec. 1981, pp 169-176, 1'auteur propose une mesuxe de 
5 densite d' echantillons par attenuation phot oni qui e en 
ligne droite au moyen de tubes de rayons X, sans 
chercher a determiner de fagon precise la dimension des 
echantillons en question. Dans les docaments 
references : Tan, S.-A. and Fwa, T". -F. , 

10 « Nondestrcutive Density Measurements of Cylindrical 
Specimens by Gamma-Ray Attenuation », JournaJ. of 
Testing Evaluation, JTEVA, Vol. 19, No. 2, March 1991, 
pp. 155-160, et : Tan, S.-A. and Fwa, T.-F., « Density 
Measurements of Cylindrical Specimens within a MoJLd by 

15 Gamma-Rays », Journal of Testing Evaluation, CTTEVA, 
Vol. 21, No. 4, July 1993, pp. 296-301, les aateurs 
proposent une mesure de densite d'echanti lions 
geologiques par attenuation photonique en ligne diroite, 
au moyen de rayonnement gamma. lis ont identifdLe et 

20 montre 1' impact des parametres geometriques des 
echantillons sur la precision de la mesure de derxsite, 
sans pour autant proposer de solution pour determiner 
avec precision lesdits parametres geometriques. 

On indique que, bien que les docixments 

25 cites ci-dessus s ' interessent a la densite des 
echantillons, il s'agit en fait de la determination de 
la masse volumique desdits objets, le terme « densite » 
etant utilise a la place de « masse volumique » par 
facilite de langage. 
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EXPOSE DE L' INVENTION 

La presente invention a pour but la 
determination de la masse volumique d' objets 
appartenant a un lot donne d' objets, par determination 
5 de la variation de la masse volumique de chacun desdits 
objets par rapport a la masse volumique, connue, de 
l'un au moins desdits objets pris comme reference ou 
etalon . 

Cette determination de la masse volumique 

10 desdits objets est effectuee au moyen d'une technique 
nucleaire non destructive, mettant en ceuvre un 
rayonnement de photons gamma, et d'une unite de 
determination par spectrometrie gamma, de 1 1 intensite 
du faisceau de photons gamma.. 

15 La determination de la masse volumique 

desdits objets passe par la determination prealable 
d'au moins une dimension significative de ceux-ci. 

Contrairement aux methodes developpees dans 
l'art anterieur cite ci-dessus, 1' invention tient 

20 compte de 1' influence des parametres geometriques , en 
mesurant de fa<?on tr§s precise au moins une dimension 
significative des objets dont on cherche a evaluer la 
masse volumique, et en utilisant cette dimension 
significative mesuree pour la determination de la masse 

25 volumique des objets controles. Ladite dimension 
significative peut §tre une largeur ou un diametre, et 
correspond a la dimension effective traversee par le 
faisceau de photons gamma. 

Le procede de determination de la dimension 

3 0 significative de l ! Qbjet fait partie du procede de 
determination de la masse volumique dudit objet. II met 
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en oeuvre une unite de mesure de dimension par 
rayonnement inf rarouge . 

On rappelle brievement que le principe 
physique de la determination de la masse volumique d'un 
5 objet par attenuation photonique consiste a irradier 
l'objet par un faisceau interrogateur de photons 
mono chromati que s d'energie E . L' intensite du faisceau 
photonique est plus ou moins attenuee en fonction de la 
masse volumique de l'objet traverse, de l'epaisseur de 
10 matiere traversee, et de la composition chimique 
elementaire de l'objet traverse, Cette intensite est 
fournie par 1' equation : 

I = I 0 exp (— p m p x) 

ou : 

15 - I est 1' intensite attenuee du faisceau 

photonique, en photons" 1 , 

- I 0 est 1' intensite non attenuee du faisceau 
photonique a l'energie E, en photons" 1 , 

- \x m est le coefficient massique 
20 d' attenuation du . photon d'energie E dans l'objet, en 



cm 2 .g-\ 



- p est la masse volumique de l'objet a 



controler, en g.cm"* 3 , 

- x est l'epaisseur de matiere traversee par 
25 le faisceau photonique, ou dimension significative de 
l'objet, en cm. 

L' expression de la masse volumique de 
1 ' objet s' en deduit directement par : 

Ln 

P = 
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Ainsi la connaissance des intensites 
transmises avec et sans interposition de 1' objet a 
controler, I et I Q respectivement, du coefficient 
massique d' attenuation yi m et de la dimension 
5 significative x de 1' objet ef f ectivement traversee 
permet de remonter a la masse volumique p dudit objet. 

La presente invention se propose de 
determiner l'epaisseur de matiere x de l'objet 
traversee par le faisceau photonique et l'intensite 

10 transmise I du faisceau photonique au niveau d'energie 
E, pour calculer ensuite la variation relative de la 
masse volumique p de cet objet par rapport a la masse 
volumique d ! au moins un objet etalon pris comme 
reference. Une caracteristique de 1' invention reside 

15 dans le fait que ces determinations d'epaisseur de 
matiere (dimension significative de 1 ! objet) et 
d' intensite du faisceau photonique sont f aites avec une 
precision micrometrique . 

La variation relative de la masse volumique 

2 0 p de 1 ! objet a controler est obtenue selon l 1 expression 
suivante : 

Ap = P - P e = 

P p • X 

ou pe est la masse volumique connue de 1 ! objet pris 
comme etalon en masse volumique, et ou x e est la 
2 5 dimension significative traversee de l'objet etalon en 
masse volumique. 

Le coefficient massique d } attenuation ]i m , 
qui depend de la composition chimique de 1 ! objet, est 



1 - 
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determine a partir d'un ou plusieurs objets etalons 
certifies et parf aitement connus, et ayant la m§me 
composition chimique que celle de I 1 objet a controler. 
II est determine lors d'une etape, qui sera decrite par 
5 la suite, d'etalonnage de 1 ' unite de determination de 
l'intensite du faisceau photonique, attenuee par la 
traversee de 1' objet etalon. 

Lorsque 1' objet a controler est de section 
circulaire, la dimension significative traversee 
10 correspond a son diametre. Lorsque 1' objet a controler 
est de forme parallelepipedique, la dimension 
significative traversee correspond a une largeur de 
1 ' ob j et . 

Par la suite on convient d'utiliser les 
15 notations suivantes lorsqu'il est necessaire de 
distinguer un objet i du lot d' objets 100 et/ou un 
objet etalon e du lot d 7 objets 100 : 

- l'indice emas est representatif des 
grandeurs relatives a un objet etalon en masse 

20 volumique, comme par exemple sa dimension significative 

x emas / 

- l'indice edim est representatif des 
grandeurs relatives a un objet etalon en dimension , 
comme par exemple sa dimension significative x e aim • 

25 Selon un premier aspect de 1 1 invention, le 

systeme pour la determination automatique de la masse 
volumique d' un objet appartenant a un lot d' objets 
comprend : 

une unite' de determination d'une 
3 0 dimension significative dudit objet, 
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une unite de determination de 
1'intensite d'un faisceau photonique, attenuee par la 
traversee dudit objet, 

une unite d' acquisition, de traitement 

5 et d 7 analyse , 

des moyens de transport de 1 ' obj et vers 
l 1 unite de determination de sa dimension significative 
et vers l 1 unite de determination de 1'intensite du 
faisceau photonique attenuee, 
lO - des premiers moyens d' ajustement de la 

position de l'objet relativement a 1' unite de 
determination de la dimension significative, 

des deux i ernes moyens d' ajustement de la 
position de l'objet relativement a 1' unite de 
15 determination de 1'intensite photonique attenuee, 

lesdits premiers et deux i ernes moyens 
d' ajustement etant capables de deplacer l'objet avec 
une precision de 1 ' ordre du micron par rapport a une 
plaque de support sur laquelle sont installes les 
2 0 elements constitutifs du systeme, 

et la position de l'objet relativement a 
1' unite de determination de 1'intensite attenuee etant 
ajustee en fonction de la dimension significative dudit 
obj et . 

2 5 De preference, 1' unite de determination 

d'une dimension significative de l'objet est une unite 
de mesure par rayonnement infrarouge. 

De preference, 1' unite de determination de 
1'intensite d'un faisceau photonique attenuee par la 

3 0 traversee de l'objet est une unite de determination par 

spectrometrie gamma, qui comprend : 
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un ensemble forme d'une source et d'un 

collimateur, 

un ensemble forme d ! un detecteur et d'un 

collimateur, 

5 - une chaine d' acquisition et de comptage 

de s photons gamma . 

L' invention met en ceuvre des moyens de 
transport, ainsi que des moyens d f ajustement de 
position de chaque objet controle par rapport a 1' unite 

10 de determination de la dimension significative de 
1' objet et/ou par rapport a 1' unite de determination de 
1' intensite attenuee, lesdits moyens d 1 ajustement de 
position etant capables de fournir une precision de 
l'ordre du micron. 

15 Selon un deuxieme aspect, 1 ' invention 

concerne un precede de mise en oeuvre du systeme pour la 
determination automatique de la masse volumique d'un 
objet (10 0) appartenant a un lot d'objets, comprend les 
etapes d'etalonnage suivantes : 

20 - une etape 1 d'etalonnage de position de 

deux ensembles infrarouges de 1 ' unite de determination 
de la dimension significative des objets, 

- une etape 2 d'etalonnage de position d'un 
support d' irradiation de 1' unite de determination par 

25 spectrometrie gamma de 1' intensite du faisceau 
photonique attenuee par la traversee des objets, 

- une etape 3 d'etalonnage de mesure d'un 
ensemble source-detecteur de 1' unite de determination 
par spectrometrie gamma de 1' intensite du faisceau 

3 0 pliotonique attenuee par la traversee des objets, 



WO 2005/100952 



PCT/FR2005/000838 



9 



et ±1 comprend des etapes de determination 
proprement dite de la dimension significative des 
objets, qui sont effectuees sur chaque objet dudit lot 
d 7 objets . 

5 Selon 1' invention, les etapes de 

determination proprement d±te comprennent : 

une etape 4 de determination de la 
dimension significative de 1' objet a controler, 

- une etape 5 de transport de 1' objet vers 
10 un support d 7 irradiation, 

une etape G d'ajustement de la position 
de 1' objet par ajustemenfc de la position du support 
d 7 irradiation par rapport a une source et un detecteur 
associe, 

15 - une etape 7 de determination de 

I'intensite attenuee du feiisceau photonique transmis a 
travers I 7 objet, 

une etape 8 d 7 acquisition, traitement et 
analyse du spectre obtenu, 

2 0 - une etape 9 de determination de la 

variation relative — de masse volumique de I 7 objet 

P 

par rapport a celle d 7 un ou plusieurs objet (s) 
etalon(s) en masse volumique, 

- une etape 10 de transport retour de 
2 5 1 7 objet jusqu' a son emplacement sur le plateau 
tournant . 

Les procedes et dispositifs selon 
l 7 invention presentent 1' interet commun d'etre rapides, 
precis, automatiques otji automat isables , et d 7 une 
30 utilisation aisee. 
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Un avantage de 1 ' invention reside dans le 
fait qu'elle associe 1 ' attenuation photonique en ligne 
droite avec une metrologie microraetrique afin de 
pallier les incertitudes liees a la meconnaissance des 
5 epaisseurs des objets traverses, et qui affectent 
directement la precision de la determination de la 
masse volumique. 

En particulier, le positionnement de chaque 
objet par rapport a l'unxte de determination de 
10 l'intensite photonique attenuee par la traversee dudit 
objet est ajuste en fonction d' une dimension 
significative dudit objet, qui a ete prealablement 
determinee par 1' unite de determination de dimension 
significative . 

15 BREVE DESCRIPTION DES DESSINS 

L 7 invention sera mieux comprise a la 
lecture de la description defcaillee qui va suivre d 1 un 
mode de realisation pref ere, donne a titre d'exemple 
non limitatif, qui est illustre par les dessins annexes 

2 0 dans lesquels : 

la figure 1 est une vue schematique de 
dessus du systeme general de determination de la 
dimension significative et de la determination de la 
masse volumique des objets, 
25 la figure 2 est une vue schematique en 

perspective du systeme general de determination de la 
dimension significative et de la determination de masse 
volumique des objets, 

les figures 3, 4, et 5 illustrent en vue 

3 0 de dessus et de fac?on schematique le dispositif de 
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determination d'une dimension significative des objets 
au moyen de rayonnement infra_rouge, et trois phases du 
procede de determination de cette dimension 
significative, 

5 la figure 6 est une vue en perspective 

et en coupe du collimateur de la source d ! irradiation 
photonique, 

la figure 7 est une vue en perspective 
et en coupe du collimateur du detecteur des photons 
10 gamma, 

la figure 8 est une vue schematique 
montrant la chaine d' acquisition et de comptage, 

les figures 9A et 9B illustrent 
1' ensemble des etapes du procede de determination de la 
15 masse volumique des objets ; la figure 9A illustre les 
etapes d'etalonnage prealables et la figure 9B illustre 
les etapes de determination proprement dites ; 

la figure 10 iLllustre la premiere etape 
du procede, qui est une etape d'etalonnage de position 
2 0 de 1 7 unite de determination de la dimension 
significative ; 

la figure 11 iLllustre la deuxieme etape 
du precede, qui est une etape d'etalonnage de position 
de 1 7 unite de determination de l 7 intensite photonique 
25 attenuee ; 

la figure 12 i Illustre la troisieme etape 
du procede, qui est une etape d'etalonnage de mesure de 
1 1 unite de determination de 1'intensite photonique 
attenuee ; 
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la figure 13 illustre la quatrieme etape 
du precede, qui est une etape de determination de la 
dimension significative d'un objet ; 

la figure 14 illustre la neuvieme etape 

5 du procede, qui est une etape de determination de la 
variation relative de masse volumique de 1' objet par 
rapport a celle d'un ou plusieurs objets etalons ; 

la figure 15 est un graphique illustrant 
la variation relative de masse volumique des objets 

D d'un lot donne d' objets par rapport a la masse 
volumique de l'un des objets pris comme etalon ou 
reference, et compare cette variation relative de masse 
volumique obtenue par 1' invention avec la variation 
relative de masse volumique theorique fournie par le 

5 fabriquant d' objets. 

EXPOSE DETAILLE D 1 UN MODE DE REALISATION PART I CUL I E R 

Les figures 1 et 2 illustrent 
respectivement en vue de dessus et en perspective un 
mode de realisation prefere du systeme general de 

0 determination par attenuation photonique de la masse 
volumique de chaque objet 100 d'un lot d' objets par 
determination de la variation relative de cette masse 
volumique par rapport a celle d' au moins un desdits 
objets pris comme reference ou etalon en masse 

5 volumique, cette determination de masse volumique 
mettant en ceuvre la determination prealable d'une 
dimension significative x dudit objet 100, et de 
1' intensity I du faisceau photonique qui irradie et 
traverse ledit objet 100. 
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Le sy st erne comprend les composants 

suivants : 

une unite 2 de determination de la 
dimension significative de l'objet 100, 
5 - une unite 3 0 de determination de 

l'intensite du faisceau photonique attenue par la 
traversee de l'objet 100, 

une unite 200 d' acquisition, de 
traitement, et d' analyse. 
10 - des moyens de transport 70, 72, 80, 82, 

84, 86, 88 et des moyens d 1 ajustement de position 74, 
76, 78, 90, 92, 94, 96, 98 de l'objet 100, 
respectivement par rapport a 1' unite 2 de determination 
de dimension et par rapport a 1' unite 3 0 de 
15 determination d'intensite, 

L ' unite 20 0 d' acquisition, de traitement et 
d' analyse est schematisee dans son ensemble a la figure 
1. Elle comprend notamment un ordinateur 170 de type PC 
sur lequel est implante un logiciel dedie qui execute 
20 des series d' instructions et des algorithmes de calcul 
du procede de determination automatique de la masse 
volumique d'objets 100 selon 1' invention. 

En se refer ant aux figures 3, 4, et 5, 
1 ' unite 2. de determination par rayonnemeiit infrarouge 
25 d'une dimension significative x de l'objet 10 0 
comprend : 

un premier ensemble infrarouge 4, 6 
const itue d'un premier emetteur infrarouge 4 et d'un 
premier recepteur infrarouge 6, 
3 0 un deuxieme ensemble infrarouge 8, 10 

constitue d'un deuxieme emetteur infrarouge 8 et d'un 
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deuxieme recepteur infrarouge 10. 

Les deux ensembles infrarouges 4, 6, et 8, 
10, sont agences de telle maniere que les axes 
respectifs 12, 14 des faisceaux de rayonnement 
5 infrarouges qu'ils engendrent sont paralleles entre 
eux, et ecartes d'une distance d. Cette distance d, 
fixee par le constructeur , est choisie de maniere a 
etre du meme ordre de grandeur que la dimension 
significative x des objets 100 a controler. Elle est 

10 ajustable. Sur l'exemple illustre, les faisceaux 
infrarouges sont orientes dans le meme sens, mais une 
configuration differente pourrait etre envisagee. 

L 7 unite 2 de determination de la dimension 
significative x de l'objet 100 par rayonnement 

15 infrarouge comprend en outre un troisieme ensemble 
const itue d'un emetteur photoelectrique 16 et d'un 
recepteur photoelectrique 18 , dispose en amont du 
premier ensemble infrarouge 4, 6, par rapport au 
deuxieme ensemble infrarouge 8, 10, le faisceau 

20 photoelectrique qu' il engendre ayant un axe 19. Sur 
1 ' exemple illustre, 1 1 axe 19 du faisceau 

photoelectrique est parallele aux axes 12, 14 des 
faisceaux infrarouges, et situe dans le meme plan que 
ceux-ci. Une configuration differente pourrait. etre 

25 envisagee, 

L 7 unite 2 de determination de la dimension 
significative x de I'objet 100 par rayonnement 
infrarouge est associee a des moyens de transport et/ou 
a des moyens d'ajustement .de position de 1'objet 100 
30 par rapport aux trois ensembles emetteurs-recepteurs 4, 
6, 8, 10, 16, 18, qui seront decrits par la suite. 
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En fonctionnement, l 1 unite 2 de 

determination de la dimension significative de l'objet 
100 se trouve dans une situation dans laquelle les 
trois ensembles emetteurs-recepteurs 4, 6, 8, 10., 16, 
5 18 sont fixes, et l'objet 100 est deplace de maniere a 
intercepter successivement le faisceau photoelectrique, 
puis le premier faisceau infrarouge, puis le deuxieme 
faisceau infrarouge . 

L ! unite 2 est gtalonnee de maniere a 

10 etablir une distance d entre les axes 12 et 14 des deux 
faisceaux de rayonnement infrarouge qui soit 
sensiblement proche de la dimension significative x e dini 
d ! un ou plusieurs objets edim pris comme etalons en 
dimension. Cet etalonnage sera decrit par la suite. II 

15 s'ensuit que lors de la determination de la dimension 
significative d f un objet 100 (non etalon) , celui-ci se 
deplace relativement aux trois ensembles emetteurs- 
recepteurs 4, 6, 8, 10, 16, 18, et passe par au moins 
une position dans laquelle il intercepte encore la 

2 0 moitie du premier faisceau infrarouge (figure 5, 
reference 22), et n' intercepte pas encore la totalite 
du deuxieme faisceau infrarouge, laissant une fraction 
du deuxieme faisceau (figure 5, reference 24) qui n'est 
pas interceptee par l'objet 100, et qui atteint le 

25 deuxieme recepteur 10 . 

La dimension significative x de l'objet 100 
se deduit de la reponse infrarouge RI correspondant a 
cette fraction de faisceau non interceptee. Cette 
dimension est obtenue par une relation du type : 

30 x = A 4 . (RI 4 ) + A 3 . (RI 3 ) + A 2 . (RI 2 ) + A x . (RI 1 ) + A 0 , 

ou A 4 , A 3/ A 2/ Ai, A 0 sont des coefficients obtenus en 



WO 2005/100952 



PCT/FR2005/000838 



16 



utilisant au mo ins quatre objets etalons en dimension 
edim, et en appliquant autant de fois cette meme 
relation, dans laquelle sont injectees la dimension 
significative x e dim connue et la reponse infrarouge RI e dim 
5 mesuree de chacun desdits objets etalons en dimension 
edim. 

Comme represents sur la figure 4, le 
troisieme ensemble 16, 18 a pour fonction de declencher 
automat iquement l'ajustement prealable de l'intensite 

10 des deux faisceaux infrarouges 22, 24, lorsqu'au cours 
du deplacement relatif selon la direction 20, l'objet 
100 intercepte le faisceau photoelectrique engendre par 
ce troisieme ensemble 16, 18. Cette operation a pour 
but de s T affranchir des perturbations environnementales 

15 telles que la salissure des optiques. Elle doit avoir 
lieu au plus tard 30 secondes avant l'operation de 
mesure proprement dite sur l'objet 100. 

La precision de la determination de la 
dimension significative x de l ! objet 100 est fonction 

2 0 de la precision du deplacement relatif de l'objet 100 
par rapport aux trois ensembles emetteurs - recept eur s 4, 
6, 8, 10, 16, 18, et done des performances et de 
1 r etalonnage des moyens de transport et/ou d 1 ajustement 
de position, aspects qui seront decrits plus en detail 

25 par la suite. 

Selon 1' invention, la determination de 
l'intensite I du faisceau photonique est effectuee au 
moyen d'une unite 3 0 de determination par spectrometrie 
gamma de l'intensite du faisceau photonique qui irradie 

30 et traverse les objets 100, comme illustre aux figures 
1, 2 et 8 . Celle-ci comprend :. 



WO 2005/100952 



PCT/FR2005/000838 



17 



un ensemble forme d'une source 
d ! irradiation photonique et d'un collimateur 32, de 
type connu en soi, 

un ensemble forme d'un detecteur et d'un 
5 collimateur 40, de type connu en soi , 

une chaine d' acquisition et de comptage 
48, de type connu en soi. 

Pour simplifier la description, la source 
d ! irradiation photonique sera par la suite nommee 
10 "source". 

Les diff<§rents elements constitutifs de 

» 

1 'unite 30 de determination sont soumis a certaines 
contraintes liees a la performance souhaitee pour le 
systeme general, ainsi qu'a l'environnement dans lequel 

15 le systeme est destine a fonctionner. Ces contraintes, 
qui concernent notamment 1' intensite de la source, le 
type de source, et les performances de la chaine 
d' acquisition et de comptage, sont les suivantes : 

- la source doit avoir une intensite telle 

2 0 que la dispersion statist ique des resultats de mesure 
soit nettement inferieure a la variation de comptage 
due a un ecart de masse volumique de l'objet a 
contr61er par rapport a la masse volumique de l'objet 
etalon pris comme reference, 

25 - l'energie de la source doit permettre un 

tres bon contraste a la suite d'une legere variation de 
la masse volumique de l'objet a controler, 

la source ne doit pas presenter une 
periode radioactive trop courte de maniere a ne pas 

30 etre contraignante en milieu industriel, 

enfin, 1' intensite et l'energie de la 
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source doivent etre compatibles avec les capacites cie 
traitement de la chaine electronique d' acquisition et 
de comptage (temps mort, empilement, saturation. . . ) . 

Selon le mode de realisation ' prefere, la 
5 source est realisee en 133 Ba d'au mo ins 10 mCi 
d'activite. Pour eviter les effets de temps mort et/ou 
de saturation, il est prefere d'utiliser une source 
dont l'activite ne depasse pas 150 mCi . La duree de 
mesure est inversement proportionnelle a l'activite de 
10 la source. 

Un exemple de realisation du collimateur 32 
de 1 7 ensemble source -collimateur, compatible avec ces 
differentes contraintes, est illustre a la figure 6. II 
comprend un coffrage de protection 34, pour assurer la 

15 securite des personnes travaillant a proximite de la 
source, qui delimite une cavite 36 dans laquelle est 
logee la source. Le faisceau de photons gamma est guide 
par une fente de collimation 38. 

Selon l 1 exemple de realisation illustre, le 

20 collimateur 32 de la source est realise en plomb et 
presente comme dimensions exterieures une hauteur de 60 
mm, une longueur de 60 mm, et une largeur de 60 mm. La 
source est une source de 133 Ba de 10 mCi d'activite, qui 
est logee dans une cavite 36 ayant un diametre de 6,1 

25 mm et une hauteur de 9,5 mm. La fente de collimation 38 
presente quant a elle une longueur de 3 0 mm, une 
largeur de 6 mm, et une hauteur de 4 mm. 

Un exemple de realisation du collimateur 40 
de 1' ensemble detecteur- col 1 imateur est illustre a la 

30 figure 7. II* comprend un coffrage de protection 42 pour 
que les rayons gamma, issus de la source et emis en 



WO 2005/100952 



PCT/FR2005/000838 



19 



dehors de la fente de collimation 38 ne soient pas 
pergus par le detecteur 49, une fente de collimation 
44, et une cavite 46 de logement du detecteur 49, 
delimitee par le coffrage de protection 42. 
5 Selon 1 1 exemple de realisation illustre, le 

collimateur 40 du detecteur 49 est realise en plomb et 
presente comme dimensions exterieures un diametre de 
14 0 mm et une longueur de 12 0 mm, et comme dimensions 
inter ieures un diametre de 80 mm et une longueur de 

10 2 00 mm. La fente de collimation presente une hauteur de 
4 mm, une largeur de 6 mm, et une longueur de 30 mm. 

Le plomb du coffrage de protection 42 peut 
etre remplace par du tungstene, qui attenue davantage 
les rayons gamma que le plomb, ce qui presente 

15 1'avantage de reduire l'epaisseur du coffrage de 
protection 42, mais le tungstene a pour inconvenient 
son prix plus el eve que celui du plomb. 

Par la suite, et pour simplifier la 
description, 1' ensemble source -collimateur sera denomme 

2 0 simplement « source » et reference 32, et 1' ensemble 

detecteur- collimateur sera denomme simplement 

« detecteur » et sera reference 40. 

La distance source -detecteur est choisie de 
fa?on appropriee. 
25 Selon le mode de realisation prefere, la 

chaine d' acquisition et de comptage 48, illustree a la 
figure 8, comprend : 

un detecteur 49, sous forme d 1 une diode 
a germanium de haute purete Ge [HP] , ayant un 

3 0 preamplif icateur , 
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un module 5 0 de traitement du signal 
(DSP : Digital Signal Processor) , 

un module haute tension 54, 
un module reseau 56 d ! acquisition et 
5 d ! interface (AIM : acquisition interface module), 

un ordinateur 170 de type PC 
d' acquisition de donnees (figure 1) . 

De fagon optionnelle la chaine 

d' acquisition et de comptage comprend un cryostat 60 
10 constitue d'un reservoir d' azote liquide, qui maintient 
a temperature constante le doigt froid de la diode Ge 
[HP], ce qui presente l'avantage de minimiser l'effet 
Doppler et d'obtenir une tres bonne resolution du 
signal, les mesures n'etant pas perturbees par un 
15 echauffement du detecteur 49. 

Le preamplif icateur est de preference 
incorpore dans la diode - Ge [HP] , ce qui presente 
l'avantage de minimiser l ! effet de capacite du au cable 
electrique et de reduire le bruit de fond electronique . 
20 II assure egalement le filtrage et la mise en forme du 
signal . 

Le signal est ensuite numerise au moyen du 
module de traitement du signal 50, puis mis en memoire. 

L' ensemble des informations obtenues 
2 5 constitue le spectre gamma, c'est a dire 1 ' histogramme 
classant dans differents canaux le nombre d' impulsions 
en fonction de leur energie . 

Les donnees sont transferees (fleche 62) 
entre le module 5 0 de traitement du signal et 
30 l 1 ordinateur 170 de l ! unite 200, d 1 acquisition, de 
traitement et d ! analyse, via le module reseau 56 
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d 1 acquisition et d ! interface, un emetteur-recepteur 63, 
et une carte reseau 59. Sur l'exemple illustre, 1 'unite 
2 00 d ! acquisition, de traitement et d ! analyse ainsi que 
la chaine 4 8 d 1 acquisition et de comptage utilisent le 
5 meme ordinateur 170, mais une configuration avec deux 
ordinateurs distincts pourrait itre envisagee. 

Cette chaine d' acquisition et de comptage 
48 est particulierement adaptee aux forts taux de 
comptage . 

10 Par ailleurs, une autre contrainte 

d' utilisation de 1' unite 30 de determination par 
spectrometrie gamma de l'intensite du faisceau 
photonique qui irradie les objets 10 0 concerne le temps 
de comptage de la chaine 48 d' acquisition et de 

15 comptage, qui doit respecter les cadences de 
fabrication des objets 100 a controler. 

Selon 1 ' invention, le temps de comptage 
peut etre une donnee d' entree du systeme, ou bien le 
resultat d' un calcul, fourni par la relation theorique 

2 0 suivante : 

2! 

A(,) 'W' p ■ r °' (rS "- i)j 

avec une approximation selon laquelle 1' angle solide 
est egal a 47CD 2 . 
ou : 

25 A(t) est 1 ' act ivite de la source, en Bq, 

D est la distance entre la source et une fenetre de 
collimation, en mm, 

s est la surface de la fenetre de collimation du 
detect eur, en mm 2 , 
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a est la largeur de l'intervalle de confiance dans 
l'hypothese ou le comptage suit une loi de Poisson, 
E est l'efficacite d' absorption totale du detecteur de 
photons , 

5 I est l'intensite du faisceau photonique & 1'energie E , 
attenuee par la traversee de l'objet, en y.s" 1 , 
I 0 est l'intensite non attenuee du faisceau photonique 
a 1'energie E, en y. s" 1 , 

Ro = — est le coefficient de transmission de l'objet 
Io 

10 traverse par les photons monochromatiques issus de la 
source , 

X est le nombre total de coups enregistres dans le 
spectre me sure , en coups, 

P est le nombre total de coups contenus dans le pic 
15 d'energie E, 

(3 sec = — est une valeur de B affectee d r un facteur de 
10 

securite egal a 10, 

avec [3 = — — , 
P 

et ou p est la masse volumique de l'objet. 

2 0 La precision de la determination de 

l f intensite I attenuee par la traversee de l T objet 100 
depend notamment de la position dudit objet 100 par 
rapport a la source 32. Elle depend done des 
performances et de 1 ! etalonnage des moyens d 1 ajustement 

25 de position. Ces aspects seront decrits plus en detail 
par la suite. 

Les differents moyens de transport 70, 72, 
80, 82, 84, 86, 88 et d T ajustement de position 74, 76, 



WO 2005/100952 



PCT/FR2005/000838 



23 

78, 90, 92, 94, 96, 98 sont illustres aux figures 1 et 
2 montrant le systeme dans son ensemble. lis ont pour 
foaction de transporter l'objet 100 vers chaque unite 
2, 3 0 de determination ou d'ajuster la position 
relative de l'objet 100 par rapport aux elements 
const itut if s de chaque unite 2, 3 0 de determination, 

Une plaque de support 15 0 supporte les 
composants du systeme general, a savoir 1' unite 2 de 
determination de la dimension significative, 1' unite 30 
de determination de l'intensite du faisceau attenue, 
les moyens de transport, les premiers moyens 
d' a. justement , et les deux i ernes moyens d' ajustement . Les 
directions de deplacement sont schematisees par le 
ref erentiel 152 de la figure 2. Les deplacements ont 
lieu dans le plan horizontal (X, Y) de la plaque de 
support 150, ou selon la direction verticale Z 
perrpendiculaire au plan horizontal (X, Y) de la plaque 
de support 150. 

Les moyens de transport 70, 72, ont pour 
0 fonction de transporter l'objet 10 0 dans une premiere 
position dans laquelle 1 1 unite 2 determine la dimension 
significative dudit objet 100. lis comprennent un 
plateau tournant horizontal 70 actionne par un moteur 
pas a pas 72, tous deux installes sur la plaque de 
5 support 150. Sur l'exemple illustre, le plateau 
tournant 70 comprend douze emplacements d'objets. 

Les premiers moyens d' ajustement 74, 76, 78 
ont pour fonction d'ajuster la. position de l'objet 100 
pair rapport aux deux ensembles infrarouges 4, 6 et 8, 
0 10, au moyen desquels est mesuree la dimension 
significative x de l'objet 100. 
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Le moyen d 1 ajustement 74 est une glissiere 
orientee suivant la direction X, dans laquelle sont 
positionnes le socle 26 de 1' unite 2 de determination 
de dimension par rayonnement infrarouge, et le plateau 
5 tournant 70. 

Les deux ensembles infrarouges 4, 6, et 8, 
10 sont installes sur le socle 2 6 de telle maniere que 
les axes 12, 14 des faisceaux infrarouges soient 
paralleles a cette direction X* Pour un lot donne 

10 d'objets dont les dimensions sont sensiblement toutes 
du meme ordre de grandeur, les positions relatives du 
socle 2 6 et du plateau tournant 70 selon cette 
direction X sont de preference f ixees une f ois pour 
toutes au debut de la serie de mesures pour le lot 

15 donne d'objets. 

Le moyen d' ajustement 76 est un verin dont 
la fonction est de rapprocher ou ecarter le premier 
ensemble infrarouge.. 4, 6 du deuxieme ensemble 
infrarouge 8, 10 suivant la direction Y. Ce deplacement 

2 0 du premier ensemble infrarouge 4, 6 selon la direction 
Y permet de positionner avec une precision 
micrometrique l'objet 100 par rapport aux deux 
faisceaux de rayonnement infrarouge pour la 
determination de sa dimension significative x, diametre 

25 ou epaisseur. 

Le moyen d' ajustement 78. est un verin dont 
la fonction est de deplacer le socle 2 6 suivant la 
direction Z. L' amplitude de ce deplacement est 
relativeTnent faible> afin d'eviter que le socle 26 

30 sorte de la glissiere 74. Le deplacement du socle selon 
la direction Z permet de connaitre avec une precision 
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micrometrique la cote de 1'objet 100 a laguelle est 
effectuee la determination de sa dimension 
significative x. 

Les tnoyens de transport 70, 72, 80, 82, 84, 
5 86, 88 ont ega lenient pour fonction de deplacer 1'objet 
100 de sa premiere position dans laquelle 1 ' unite 2 
determine la dimension significative x vers sa deuxieme 
position dans laquelle l 1 unite 30 determine l'intensite 
attenuee I du faisceau photonique. lis comprennent le 

10 plateau tournant 70 actionne par son moteur pas a pas 
72. En effet, plusieurs objets 100 etant disposes en 
cercle sur le plateau tournant 70, la rotation dudit 
plateau 70 remplit deux actions simultanees consistant 
d 7 une part a transporter un objet 100 vers sa premiere 

15 position de mesure et d 7 autre part a eloigner 1' objet 
100 precedent de sa premiere position de mesure pour 
1 7 emmener vers une position intermediaire apres lui 
avoir fait effectuer un deplacement angulaire d 7 angle 
A. Sur 1 ' exemple illustre aux figures 1 et 2 , cet angle 

20 A est de 90°. Les moyens' de transport comprennent 
egalement un bras de manutention 80 qui saisit l 7 objet 
100 installe suit le plateau tournant 70 en sa position 
intermediaire et le transporte sur un support 
d 7 irradiation 90 dispose entre le collimateur 32 de la 

25 source et le collimateur 4 0 du detecteur. Selon 
1 7 exemple illustre a la figure 2, le bras de 
manutention 80' comprend une pince de prehension 82 
articulee sur urn segment intermediaire 84, lui-meme 
articule sur un verin 86 capable de se deplacer en 

3 0 translation selon la direction X de la plaque de 
support 150, le long de rails de guidage 88. Les 
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raouvements de serrage/desserrage de la pince 82 et de 
pivotement de celle-ci par rapport au segment 84, ainsi 
que les mouvements de pivotement du segment 84 par 
rapport au verin 86 sont commandes par des actionneurs 
5 (non representes) . 

Les deux i ernes moyens d'ajustement 90, 92, 
94, 96, 98 ont pour fonction d'ajuster la position de 
1' objet 100 par rapport a. la source 32 et au detecteur 
4 0 de 1' unite 3 0 de determination par spectrometrie 

10 gamma de l'intensite du faisceau qui va traverser ledit 
objet 100. lis comprennexit le support d' irradiation 90 
sur lequel est installe 1 ' objet 10 0. Ce support 
d' irradiation 90 presente une face superieure 92 a 
section en forme de V, ou tout autre moyen equivalent 

15 pour que 1' objet 100 se fcrouve automat iquement installe 
en position d'equilibre stable sur ledit support 
d ! irradiation 90, et notamment qu'il soit empeche de se 
deplacer par rapport au support d' irradiation 90 selon 
la direction X de la plaque de support 150. Le support 

20 d' irradiation 90 est positionne suivant la direction X 
de la plaque de support 150 au moyen d' une glissiere 
94, qui, de preference, est confondue ayec la glissiere 
74. Pour un lot donne d'objets 100, ce positionnement 
est effectue une fois pour toutes au debut de la serie 

25 de mesures correspondant a un lot donne d f objets. Le 
support d' irradiation 90 peut etre deplace selon la 
direction Y de la plaque de support 150 au moyen d'un 
verin 96 et selon la diirection Z perpendiculaire a la 
plaque de support 150 au moyen d'un verin 98. Les 

3 0 ajustements de position effectues par les verins 96 et 
98 permettent de centrer* sensiblement 1' objet (suivant 
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la direction Z) entre les fentes des collimateurs 
respectifs de la source 32 et clu detecteur 40. 

D' autre part il est necessaire de 
positionner avec une precision micrometrique l'objet 
5 suivant la direction Y de telle maniere que la mesure 
d' intensite 1 du faisceau photonique s'effectue 
exactement a la cote de l'objet a laquelle a ete 
effectuee la determination de sa dimension 
significative x. Ce positionnement est effectue par une 

10 mise en butee de l'objet sur la face superieure 92 du 
support d' irradiation 90. Cette mise en butee peut etre 
effectuee par exemple par une operation de soufflage au 
raoyen d'un dispositif de soufflage (non represente, qui 
envoie de l'air comprime star l'objet, suivant la 

15 direction Y, afin de plaquer celui-ci contre une butee 
93 du support d' irradiation 90, 

Sur la figure 1. sont illustrees les 
liaisons par des moyens de connexion appropries 180 
entre d'une part les different s verins 76, 78, 86, 96, 

2 0 ' 98 de deplacement des pieces mobiles en translation et 

le moteur pas a pas 72 qui entraine en rotation le 
plateau tournant 70, et d' autre part des unites 160 de 
control e et de pilotage. Ces unites 160 commandent la 
me c an i que et 1 ' automatique du systeme, et sont reliees 
25 a l 1 unite centrale 172 de 1 ' ordinateur 170 de 1' unite 
200 d' acquisition, de traitement et d' analyse, par 
d'autres moyens de connexion appropries 190. 

On va maintenant decrire le procede de 
determination de la masse volumique p de chacun des 

3 0 objets d'un lot donne d'objets 100, par comparaison 

avec la masse volumique p e mas d'un ou plusieurs objet(s) 
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choisi(s) comme etalon(s) ou reference (s) en masse 

volumique, et faisant partie du meme lot d' objets 100. 

Le procede est mis en ceuvre au moyen 

d'algorithmes traduisant des series d # instructions qui 
5 perraettent d'effectuer de maniere automatique les 

dif ferentes etapes du procede. 

Le procede de 1 ' invention comprend des 

etapes prealables d' etalonnage, qui sont effectuees une 

fo±s pour toutes avant de commence r une serie de 
10 itiesures sur un lot donne d' objets, et des etapes de 

determination proprement dites, qui sont effectuees sur 

chaque objet 100 dudit lot d'objets- L'ensemble des 

etapes du procede est represents de f agon synthetique 

sur les figures 9A et 9B . 
15 Les etapes d' etalonnage respectent une 

chronologie pre-etablie, et concernent les composants 

suivants du systeme : 

- etape 1 : etalonnage de position des deux 

ensembles infrarouges 4, 6 et 8, 10 de 1 ' unite 2 de 
20 determination de la dimension significative des objets 

10O, 

etape 2 : etalonnage de position du support 
d' irradiation 90 de 1' unite 30 de determination par 
spectroraetrie gamma de l'intensite du faisceau 
25 photonique attenuee par la traversee des objets 100, 

etape 3 : etalonnage de mesure de 
l'ensemble source-detecteur 32, 40 de 1' unite 30. 

L' etape 1 d'etalonnage de position des deux 
ensembles infrarouges 4, 6 et 8, 10 est illustree sur 
3 0 la figure 10. 
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Cette etape 1 d'etalonnage consiste a 
ajuster la position suivant la direction Y du premier 
ensemble infrarouge 4, 6 par rapport au deux i erne 
ensemble infrarouge 8, 10, afin de fixer 1 7 ecartement d 
5 entre les faisceaux infrarouges emis respect ivement par 
Xes deux emetteurs 4, 8, en fonction de la dimension 
significative x e dim connue precisement d'un ou plusieurs 
objet (s) etalon(s) en dimension edim. En pratique, 
1' ecartement d est etabli en eloignant progressivement 

10 le premier ensemble infrarouge 4, 6 du deuxieme 
ensemble infrarouge 8, 10 suivant la direction Y, ce 
dernier restant fixe a une position Y F i X , et en mesurant 
l_a reponse infrarouge de 1' objet pour c ha que position 
clu premier ensemble infrarouge 4, 6. 

15 L' etape 1 d'etalonnage de position des deux 

ensembles infrarouges 4, 6 et 8, 10 comprend tout 
dL'abord une saisie par 1 ' operateur d'un ensemble de 
parametres d' entree au moyen d'un module interact if . 
Ces parametres comprennent : 

2 0 - la configuration des coraposants qui ont 

um deplacement micrometrique : verins 76, 78, ce qui 
permet de gerer leur cinematique : position, vitesse, 
acceleration, 

- la configuration du plateau tournant 70, 
25 o ' est-a-dire la nature des objets qui occupent les 
cliff erents emplacements sur le plateau tournant 7 0 : 
objet quelconque 10 0, ou objet etalon en dimension 
edim, ou objet etalon- en masse volumique emas, ou 
emplacement libre, 

3 0 - 1 ' emplacement occupe par les objets 

etalons en dimension edim sur le plateau tournant 70, 
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cet emplacement etaxit un nombre allant de 1 a 12 pour 
1 ' exemple illustre , 

- la position Z meS ure suivant la direction Z 
du socle 26 de l'unite 2, qui correspond a une cote 
5 Zedim sur 1 ' objet edim par rapport a la base de celixi- 
ci, 

les positions Y (1) et Y (N) bornant 
l'intervalle de deplacement du premier ensemble 
infrarouge 4, 6 suivant la direction Y, 
10 - le pas INT exprime en \\m du deplacement 

du premier ensemble infrarouge 4, 6 suivant la 

Y - Y 

direction Y, (il f aut que DEP soit un nombre 

INT 

entier) . 

L'etape 1 d'etalonnage de position des deux 
15 ensembles infrarouges 4, 6 et 8, 10 comprend ensuite 
les operations automatisees suivantes : 

a) dep3_acement du - - socle 26 suivant la 
direction Z jusqu'a la position Z meS ure/ par actionnement 
du verin 78, 

2 0 b) deplacement angulaire du plateau 

tournant 7 0 afin de transporter 1' objet etalon en 
dimension edim juscjuL'a sa position initiale de mesure 
par rapport a 1' unite 2, 

c) deplacement du premier ensemble 
25 infrarouge 4, 6 siaivant la direction Y jusqu'a sa 

position de depart Yd), par actionnement du verin 76, 

d) deplacement progressif, par increments 
successifs de INT, du premier ensemble infrarouge 4, 6 
suivant la direction Y, en l'eloignant du deux i feme 

3 0 ensemble infrarouge 8, 10 fixe a une position Y F:r x/ 
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entre les positions Y(l) et Y (N) , et determination 
simultanee de la reponse infrarouge RI (n) de l'objet 
edim, correspondant a chaque position Y(n) , de la 
maniere suivante : 
5 &-1) depl a cement angulaire du 

plateau tournant 7 0 afin de transporter l'objet etalon 
en dimension edim ^. sa position finale de mesure, 

dL-2 ) mesure de la reponse infrarouge 
RI (n) dudit objet etalon en dimension edim, 
10 <i-3) deplacement angulaire du 

plateau tournant 70 afin de ramener l'objet etalon en 
dimension edim a sa position initiale de mesure, 

e) calcul de la reponse infrarouge optimale 

1X1 opt — 2 

15 ou : RImin est la valeur de la saturation minimale de la 
reponse infrarrouge : au debut de 1 ' etalonnage , 
l'ecartement des deux ensembles infrarouges 4 7 6, 
et 8, 10 est tres inferieur a la dimension 
significative x e dim de l'objet etalon en dimension 

2 0 edim ; par suite, lorsque 50% du premier faisceau 

infrarouge est intercepts par l'objet edim, 100% 
du deuxieme faisceau infrarouge est intercepts par 
cet objet edim ; les premieres reponses 
infrarouges oat alors une valeur identique RImin 

25 dite « saturSe », 

et : RImax est la vaZLeur de la saturation maximale de la 
reponse infrarouge : a la fin de 1 ' etalonnage, 
l'ecartement des deux ensembles infrarouges 4, 6, 
et 8, 10 est tres superieur a la dimension 

30 significative :x e dim de l'objet etalon en dimension 

edim ; par suite, lorsque 50% du premier faisceau 
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infrarouge est intercepts par l'objet edim, 0% du 
deuxieme faisceau infrarouge est intercepts par 
cet objet edim ; les dernieres reponses 
infrarouges ont alors une valeur identique RImax 
5 dite « saturee », 

f) calcul de la position optimale Y 0 pt du 
premier ensemble infrarouge 4, 6 par rapport au 
deuxieme ensemble infrarouge 8, 10 : la reponse 
infrarouge optimale RIopt est comprise entre deux 
10 valeurs successives RI(j) et RI (k) precedemment 
calculees de la reponse infrarouge, qui correspondent 
respect ivement a deux positions Y(j) et Y(k) du premier 
ensemble infrarouge 4, 6 ; la position optimale YOPT 
s'en deduit de la maniere suivainte : 
RI OPT -RI(j) 

15 S± T) T T,T/tN <1 ' al0rS Y OPT = 

RIopt -Ri(k) 
RI opr -RI(j) 

si >1 , alors Y 0 pt = Y(k) 

RI OPT -RI(k) 

Les operations a) a f) ci-dessus peuvent 
etre renouvelees avec autant d r autres objets etalons en 
dimension edim que necessaire. 

2 0 A 1 ' issue de l'etape 1 d'etalonnage de 

position des deux ensembles infrarouges 4, 6 et 8 7 10, 
il est cree un premier fichier d'etalonnage qui 
comporte, notamment, 1 ' ecartement optimal d des deux 
ensembles infrarouges 4, 6 et 8, 10 correspondant 

25 sensiblement a la dimension significative des objets 
d = |[y F ix - Y 0 pt|| • 

L'etape 2 d'etalonnage de position du 
support d' irradiation 90 de 1' unite 30 de determination 
par spectrometrie gamma de IL'intensite du faisceau 
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photonique attenuee par la traversee d^s objets 10 0 est 
illustree de fagon synthetique sur la figure 11. 

Cette etape 2 d' etalonraage consiste a 
ajuster la position suivant la direction Z du support 
5 d' irradiation 90 par rapport a la source 32 et au 
detecteur 40 associe, afin de fixer la position Zopt 
suivant la direction Z de la face superieure 92 du 
support d' irradiation 90 sur lequel sont positionnes 
les objets 100 traverses par le faisceau photonique, en 

10 f one t ion de la masse volumique p connue precisement 
d'un ou plusieurs objet(s) etalon (s) en masse volumique 
emas. En pratique, la position Zqpt est etablie en 
deplagant progressivement le support d' irradiation 90 
suivant la direction Z, et en irradiamt plusieurs fois 

15 l'objet etalon en masse volumique emas installe sur le 
support d' irradiation 90 pour chacune des positions de 
celui-ci. Elle resulte du calcul d' am minimum d'une 
regression polynomiale d'ordre 4. EILle comport e une 
etape de determination de la dimension significative 

2 0 Xemas de chaque objet emas etalon en masse volumique. 

L 7 etape 2 d' etalonnage de position du 
support d' irradiation 90 de 1' unite 3 0 comprend tout 
d'abord une saisie par l'operateur d'un ensemble de 
parametres d 7 entree au moyen d'un module interactif . 

2 5 Ces parametres comprennent : 

- la configuration des compos ants qui ont 
un deplacement micrometrique : verins 96, 98, ce qui 
permet de gerer leur cinematique : position, vitesse, 
acceleration, 

30 - la configuration du - plateau tournant 70, 

e'est-a-dire la nature des objets <qui occupent les 
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differents emplacements sur le plateau, tournant 70 : 

objet quelconque 10 0, ou objet etalon en dimension 

edim, ou objet etalon en masse volamique emas, ou 

emplacement libre, ■ < 

5 - 1' emplacement occupe par les objets 

etalons en masse volumique emas sur le plateau tournant 

70, cet emplacement etant un nombre alLlant de 1 a 12 

pour 1'exemple illustre, 

- la duree de mesure ou temps de comptage, 
10 - les positions Z(l) et Z (N) bornant 

l'intervalle de deplacement du support d' irradiation 90 

suivant la direction Z, 

le nombre M de mesures de 1 ' intensite 

photonique attenuee par la traversee c3.e 1' objet, pour 
15 chaque position Z(i) occupee pax* le support 

d' irradiation, pour i = 1,~.,N. 

L'etape 2 d'etalonnage c3.e position du 

support d' irradiation 90 de 1 ' unite 30 comprend ensuite 

les operations automatisees suivantes : 
20 a) determination de la dimension 

significative x em as de 1 ' objet etalon en masse 

volumique, conformement a l'etape 4 cjui sera decrite 

ci-apres , 

b) deplacement angulaire du plateau 
25 tournant 70, d'un angle A, afin de transporter 1' objet 

etalon en masse volumique emas en une position 
intermediaire ou il est saisi par le biras de prehension 
80, 

c) positionnement de 1' objet emas sur le 
30 support d' irradiation 90, qui comporte les sous- 
operations suivantes : 
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c-1) deplacement du support 

d' irradiation 90, suivant la direction Z et venrs le 
has, par actionnement du verin 98, 

c-2) deplacement du bras de 
5 manutention 80 depuis sa position d'attente jusqix'a la 
verticale de la position intermediaire de l'objet emas, 
par actionnement du verin 86, 

c-3) saisie, par le bras de 
manutention 80, de l'objet emas, et transport de celui- 
10 ci jusqu'a la verticale de la face superieure 92 du 
support d' irradiation 90, par actionnement du verin 86, 

c-4) deplacement du support 

d' irradiation 90 jusqu'a la position Z(l), suiva_nt la 
direction Z et vers le haut, par actionnement du verin 
15 98, 

c~5) depose de l'objet emas s~ur la 
face superieure 92 du support d' irradiation 90, par le 
bras de manutention 80, par actionnement du verin S6,. 

c-6) deplacement retour du bras de 
20 manutention 80 jusqu'a sa position d'attente,. par 
actionnement du verin 86, 

c-7) mise en butee de l'objet emas 
sur la face superieure 92 suivant la direction Y~ r par 
exemple par une operation de soufflage, qui se deroule 
25 de la maniere suivante : 

deplacement du • support 
d' irradiation 90 suivant la direction Z et vers l_e bas 
jusqu'a une position dite de soufflage, dans lacguelle 
l'objet se trouve en regard d'un dispositi f de 
3 0 soufflage prevu dans le systeme, 
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envoi d'air comprime 
provenant du dispositif de soufflage sur l'objet emas r 
suivant la direction Y, afin de plaquer celui-ci contre 
une butee 93 du support d 7 irradiation 90, 
5 d) ajustement proprement dit de la position, 

du support d 7 irradiation 90 par rapport a la source 32 
et au detecteur associe 40, qui comporte les sous — 
operations suivantes : 

d-1) deplacement progress if du support 
10 d 7 irradiation 90 suivant la direction Z entre la 
position Z(l) predeterminee et la position Z (N> 
predeterminee , 

d-2 ) pour chaque position Z(i), 
i = 1, N, irradiation de 1'objet etalon en masse 

15 volumique emas par le faisceau photonique, un nombre M 
de fois, ce qui conduit a un ensemble de valeurs 
d'intensite attenuee I(i, j ) , ou i = 1, N represents 
le nombre de positions successives Z(i) occupees par le 
support d' irradiation 90 et j = 1, M represente le 

20 nombre d 7 irradiations effectuees a chaque position 
Z(i), 

d-3) calcul de la position optimale 
Zqpt du support d 7 irradiation 90 a partir d'une 
regression polynomiale d'ordre 4 des positions Z(i) pax 
25 rapport aux intensites attenuees I(i, j ) , cette 
regression polynomiale d'ordre 4 etant predeterminee et 
integree comme une donnee de l 7 unite 200 d' acquisition j 
de traitement et d 7 analyse, 

e) transport retour de l 7 objet etalon e:n 
3 0 masse volumique emas sur le plateau tournant 70, pazr 
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une sequence d/ operations inverses aux sous -operations 
c-1) a c-6) qui ont ete detaillees ci-dessus. 

A 1 ' issue de l'etape 2 d'etalonnage de la 
position du support d' irradiation 90 de l'unite 30 de 
5 determination par spectrometrie gamma de 1'intensite du 
faisceau photonique attenuee par la traversee des 
objets 100 il est cree un deuxieme fichier d'etalonnage 
qui comporte, notamment, la position optimale Zqpt du 
support d' irradiation 90 suivant la direction Z. 
10 L'etape 3 d'etalonnage de mesure de 1' unite 

3 0 de determination par spectrometrie gamme comprend 
les operations automatisees suivantes : 

a) determination de 1'intensite photonique 
I emas attenuee par la traversee d' un objet etalon en 

15 masse volumique emas pris comme reference, 

b) calcul du coefficient massique 
d' attenuation JI m de I 7 objet etalon en masse volumique, 
et par suite de tous les objets du lot d' objets par la 
relation : 

on n — - T ^ emas 

20 Pemas " n T * 

r 4, m emas 0 

A 1' issue de l'etape 3 d'etalonnage de 
mesure de 1' unite 3 0 de determination par spectrometrie 
gamma de 1'intensite du faisceau photonique attenuee 
par la traversee des objets 100 7 il est cree un 

25 troisieme fichier d'etalonnage qui comporte, notamment, 
1'intensite photonique I e mas attenuee par la traversee 
de 1' objet etalon en masse volumique emas. 

Les etapes de determination proprement 
dites respectent elles aussi une chronologie pr<§- 

3 0 etablie et concernent les operations suivantes : 
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etape 4 : determination de la dimension 
significative x de 1'objet 100 a controler, 

etape 5 : transport de 1'objet 100 vers 
le support d' irradiation 90, 
5 - etape 6 : ajustement de la position de 

1'objet 100 par ajustement de la position du support 
d' irradiation 90 par rapport a la source 32 et au 
detecteur associe 40, 

etape 7 : determination de l'intensite 
0 attenuee I du faisceau photonique transmis a travers 
1'objet 100, 

etape 8 : acquisition, traitement et 
analyse du spectre obtenu, 

etape 9 : determination de la variation 

5 relative de masse volumique de 1'objet 100 par 

P 

rapport a celle d'un ou plusieurs objet(s) etalon(s) en 
masse volumique emas, 

etape 10 : transport retour de 1'objet 
100 jusqu'a son emplacement sur le plateau tournant 70. 

0 L' etape 4 de determination de la dimension 

significative x de 1'objet 10 0 a controler est 
illustree de fa?on synthetique a la figure 13 . Elle 
comprend tout d'abord une saisie par l'operateur d'un 
ensemble de parametres d' entree au moyen d'un module 

5 interactif . Ces parametres comprennent : 

- la configuration des compos ants qui ont 
un deplacement micrometrique : verins 76, 78, ce qui 
permet de gerer leur cinematique : position, vitesse, 
acceleration, 
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- la configuration du plateau tournant 70, 
c'est-a-dire la nature des objets qui occupent les 
different s emplacements sur le plateau tournant 7 0 : 
objet quelconque 10 0, ou objet etalon en dimension 

5 edim, ou objet etalon en masse volumique etnas, ou 
emplacement libre, 

- 1 7 emplacement occupe par 1' objet 100 sur 
le plateau tournant 70, cet emplacement etant un nombre 
allant de 1 a 12 pour l'exemple illustre, 

10 - la position Z me sure suivant la direction Z 

du socle 26 de 1' unite 2, qui correspond a une cote z 
sur 1 7 objet 10 0 par rapport i. la base de celui-ci, 

- le nombre P de mesures infrarouges pour 
chaque objet etalon en dimension edim(n) 7 n = 1, N, 

15 ou N est le nombre d 7 objets etalons en dimension, 

- le nombre Q de mesures infrarouges pour 

1' objet 100. 

L'etape 4 de determination de la dimension 
significative x de 1 7 objet 100 a contrdler utilise 
20 aussi les donnees contenues dans le premier fichier 
d'etalonnage issu de l'etape 1. 

L'etape 4 de determination de la dimension 
significative x de 1 7 objet 100 a controler comprend 
ensuite les operations automatisees suivantes : 
25 a) deplacement du socle 2 6 suivant la 

direction Z jusqu'a la position Zmesure, pzar 
actionnement du verin 78, 

b) deplacement du premier ensemble 
infrarouge 4, 6 suivant la direction Y, par 
3 0 actionnement du verin 76, jusqu 7 a une position Y mesu re 
def inie par : 
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Y m esure = Yqpt + (x e dim ~~ X e dimMOY) 

OU : 

Yopt est la position optimale provenant de 
1'etape 1 d' etalonnage, cette valeur etant contenue . 
5 dans le premier fichier d' etalonnage , 

x e di m est la dimension de l'objet etalon en 
dimension edim utilise lors de l'Stape 1 d' etalonnage, 
cette valeur etant contenue dans le premier fichier 
d' etalonnage, 

10 XedimMOY est la dimension significative 

moyenne de tous les objets edim e talons en dimension, 
cette valeur etant fournie par le fabricant, 

c) mesure repetee P fois de la reponse 
infrarouge RI (p) , p = 1, P des N objets etalons en 

15 dimension edim(n) , n = 1, ... N, ce qui conduit a un 
ensemble de valeurs RI (n, p) , 

d) calcul proprement dit de la dimension 
significative x de l'objet 100 de la maniere suivante : 

VRI(n } p) 

d-1) calcul de la moyenne RI edimMOY = — — 

2 0 des reponses infrarouges de chaque objet etalon en 
dimension edim(n) dont la dimension significative 
Xedim (n) est connue, et mise en oeuvre d'une regression 
polynomiale d'ordre 4 des dimensions signif icatives 
Xedim (n) pour calculer les coefficients A 0 , Ai, A 2 , A 3 , A 4 

25 d 7 une relation du type : 

Xedim (n) = A 4 . (RIedimMOY (n) ) 4 + A 3 . (RIedimMOY (n) ) 3 + 

A 2 . (RIedimMOY (n) ) 2 + Ai . (RIedimMOY (n) ) 1 + A 0 , 

d-2 ) mesure, repetee Q fois, de la reponse 
infrarouge RI (q) , q = 1, Q de l'objet a controler 
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10 0, calcul de la moyenne RI = — de ces reponses 

infrarouges, et calcul de la dimension significative x 
recherchee de 1'objet 100 par la relation : 
x = A 4 . (RI) 4 + A 3 . (RI) 3 + A 2 . (RI) 2 + Ai . (RI) 1 + A 0 
5 L'etape 5 de transport de 1'objet 100 a 

controler vers le support d' irradiation est une etape 
automat isee qui reprend la sequence des sous -operations 
b) et c) de 1' etape 2 d'etalonnage qui ont ete 
detaillees ci-dessus . 

10 L' etape 6 d'ajustement de la position de 

1'objet 100 par rapport a la source 32 et au detecteur 
associe 4 0 est une etape automat isee qui reprend la 
sous -operation d) de 1' etape 2 d'etalonnage qui a ete 
detaillee ci-dessus . 

15 L' etape 7 de determination de 1' intensity 

photonique I du faisceau photonique, attenuee par la 
traversee de l'objet 100 consiste en une mesure 
d'activite, qui est ensuite acquise, traitee et 
interpretee de maniere connue en soi • 

20 L ' etape 8 d' acquisition, de traitement et 

d' analyse du spectre obtenu est une etape automatisee 
qui met en oeuvre des algorithmes de calculs connus en 
soi qui sont executes par le logiciel dedie imp ante sur 
1'ordinateur 170 de 1' unite 2 00 d' acquisition, de 

25 traitement et d 7 analyse. 

L' etape 9 de determination de la variation 

relative de masse volumique de l'objet 100 par 

P 

rapport a celle d'un ou plusieurs objet(s) etalon(s) en 
• masse volumique emas est illustree de fagon synthetique 
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a la figure 14. C'est une etape automatisee de calcul 
dans laquelle interviennent 1' equation 

I _ 

et les donnees qui ont £t<§ 



P 



a I 

1 ± 



M* m P emas ^ emas 



deterttilnees aux differentes etapes precedentes . 
5 L' etape 10 de transport retour de 1'objet 

10 0 suit son emplacement du plateau tournant 7 0 est une 
etape automatisee qui rep rend la sous -operation e) de 
1 ' etape 2 d'etalonnage qui a ete detaillee ci-dessus. 

Le procede qui vient d'etre decrit est mis 

10 en ceuvre au moyen d'un logiciel dedie. Ce logiciel 
comporte cinq modules independants, et un menu 
principal interactif, au moyen duquel un operateur 
choisit de faire executer l'un des cinq modules. Les 
cinq modules regroupent les fonctions suivantes : 

15 - premier module : determination de la 

masse volumique d'un objet, qui recouvre 1' etape 3 
d'etalonnage et les etapes 4 a 10 de determination 
proprement dite de la masse volumique, 

deuxieme module : determination de la 

2 0 dimension significative d'un objet, 

- troisieme module : etalonnage de position 
de 1' unite de determination de dimension significative, 

- quatrieme module : etalonnage de position 
de 1 'unite de determination de 1' intensity photonique 

25 attenu.ee, 

cinquieme module : gestion des fichiers 

de doixnees . * 
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EXEMPLE 

Le systdme et le procede decrits 
precedemtnent ont ete testes. 
5 La source etait une source de 133 Ba de 10 

mCi d'activite. La duree des acquisitions etait de 
l'ordre de 20 minutes. 

Les mesures ont ete effectuees sur un lot 
de 7 pastilles d'oxyde d' uranium (U0 2 ) dont les 
10 caracteristiques : diametre, hauteur, masse volumique 
ou densite, apparaissent dans le tableau I : 



N° pastille : i 


1 


2 


3 (etalon) 


Diametre (mm) 


8, 165 


8,143 


8, 166 


Hauteur (mm) 


11, 54 


11,44 


11, 27 


Masse volumique 
(g.cm 3 ) ; 


10,260 
± 0.003 


10, 130 
± 0.003 


9, 900 
± 0.003 


ecart-type 
(g . cm 3 ) 


1, 99 . 10" 2 


1,98.10" 2 


1,96.10" 2 



Tableau I 



N° pastille : i 


4 


5 


6 


7 


Diametre (mm) 


8, 147 


8, 123 


8,117 


8, 169 


Hauteur (mm) 


11,49 


11, 29 


11, 54 


11, 59 


Masse 
volumique 
(g.cm 3 ) 


10, 150 
± 0.003 


9, 950 
± 0.003 


9, 960 
+ 0.003 


10, 070 
± 0.003 


ecart-type 
(g.cm 3 ) 


1,98.10" 2 


1,95.10' 2 


1, 93 . 10" 2 


1, 97.10" 2 



15 



Tableau I (suite) 
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La pa_stille 3 est prise comme pastille 

etalon. 

L'objectif des mesures est la determination 
precise, au moyen du systeme et du precede de 
1' invention, de la variation relative de masse 
volumique des pastilles (1,2,4,5,6 et 7) par rapport a 
la masse volumique de la pastille etalon (pastille 3) . 
La relation suivante s' applique : 



Ap 
P 



x 

x. 



1- 
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10 Les diametres des pastilles supposees 

«inconnues » sont obtenus par l'etape de determination 
de la dimension significative, ici le diametre . des 
pastilles, par rayonnement inf rarouge . 

Les resultats des comptages obtenus par 

15 spectrometrie gamma, pour chacune des six pastilles 
sont exposes dans le tableau II. lis ont ete obtenus en 
suivant scrupuleusement la chronologie du procede 
presentee precedemment . 



N° 

PASTILLE 


I (en coups) 


E CARTS DE MASSE VOLUMIQUE 


1 


974725 ± 1974 


(3,448633 ± 0,017045). 10" 2 


2 


1012550 ± 2012 


(2,286541 ± 0,061460). 10~ 2 


4 


1009661 ± 2010 


(2,344449 ± 0,016572). 10" 2 


5 


1063886 ± 2063 


(6,611105 ± 0,132441). 10~ 3 


6 


1067853 ± 2067 


(5,941459 + 0,122442). 10" 3 


7 


1014895 ± 2015 


(1,873675 ± 0,017101). 10" 3 
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Tableau II 

Les ecarts types relatifs aux variations de 
masses volumiques mesurees ont ete estimes par un 
5 calcul de propagation d' incertitudes . Le tableau III 
est un tableau compa.ratif de ces resultats avec les 
ecarts theoriques fournis par le fabriquant de 
pastilles . 



N° Pastille 


1 


2 


4 


[Ap/p]theorique 


3, 63636 . 10" 2 
± 2,77. 10' 3 


2, 32323 . 10" 2 
± 2,78.10" 3 


2, 52525. 10" 2 
± 2, 77. 10" 3 


[Ap/p] meS ur£ 


3,44863 . 10" 2 
± 1,70.10 -4 


2, 28654 . 10" 2 
± 1,65.10" 4 


2, 34445. 10" 2 
± 1,66.10" 4 


Tableau III 


N° Pastille 


5 


6 


7 


[Ap/p]theori 
que 


5, 0505. 10 -3 
± 2, 79 . 10" 4 


6,0606.10~ 3 
± 2,77.10" 4 


1, 71717. 10" 2 
± 2,76.10" 3 


[Ap/p]mesure 


6, 611. 10" 3 
± 1, 32 .10~ 4 


5, 9415. 10" 3 
± 1,22.10" 4 


1, 8737 . 10" 2 
+ 1,71.10" 4 



Tableau III (suite) 

15 

Ces resmltats sont illustres par le 
graphique de la figure 15 . Les cercles representent les 

valeurs de resultants de la mesure, alors que les 

P 
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croix representent les valeurs cie fournis par le 

P 

fabricant. L' intervalle materialise represente l'ecart- 
type calcule a partir cie donnees fournies par le 
f abricant • 

5 Ces resultats montrent que le systeme et le 

procede selon 1' invention permettent de detecter une 
variation relative de masse volumique d' environ 6.10" 3 
par rapport a la pastille choisie comitie objet etalon. 
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RE VEND I CAT I ONS 

1. Systeme, pour La determination 

automatique de la masse volumique d'un objet (100) 
5 appartenant a un lot d'objets, caracterise en ce qu'il 
comprend : 

une unite (2) de determination d'une 
dimension significative (x) dudit objet (100), 

une unite (3 0) de determination de 
10 l'intensite (I) d'un faisceau photonique, attenuee par 
la traversee dudit objet (10 0) , 

une unite d 1 acquisition, de traitement 
et d 1 analyse (2 00), 

des moyens de transport (70, 72, 80, 82, 
15 84, 86, 88) de l'objet (100) vers l'unite (2) de 
determination de sa dimension signif icative (x) et vers 
1' unite (3 0) de determination de l'intensite photonique 
attenuee , 

des premiers moyens d'ajustement (74, 
20 76, 78) de la position de 1' objet (100) relativement a 
1' unite de determination de sa dimension significative 
(x) , et 

des deuxiemes moyens d'ajustement (90, 
92, 94, 96, 98) de la position de l'objet (100) 
25 relativement a 1' unite (30) de determination de 
l'intensite photonique attenuee, 

en ce que lesdits premiers moyens 
d'ajustement et deuxiemes moyens d'ajustement sont 
capables de deplacer l'objet (100) avec une precision 
30 de 1'ordre du micron par rapport a une plaque de 
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support (150) sur laquelle sont installes les elements 
const itut if s du systeme, 

et en ce que la position de 1' objet (100) 
relativement a 1 ' unite (3 0) de determination de 
5 1' iLntensite attenuee (I) est ajustee en fonction de la 
dimension significative (x) dudit objet (100) . 

2. Systeme selon la revendication 1, 
caracterise en ce que 1 'unite (200) d 1 acquisition, de 
traitement et d 1 analyse comprend un ordinateur (170) 

10 dans lequel est imp 1 ante un logiciel ded±e qui execute 
des series d' instructions et des algorithmes de calcul 
automatique de la masse voluraique de l'objet (100). 

3. Systeme selon la revendication 1 ou 2, 
caracterise en ce que l 1 unite (200) d ? acquisition, de 

15 traitement et d 1 analyse fournit la variation relative 

(^^) de la masse volumique (p) de l'objet (100) par 
P 

rapport a celle, connue, d'au moins un objet etalon 
(emas) appartenant au meme lot d'objets (100) . 

4. Systeme selon 1 ' une quelconque des 
20 revendications 1 a 3, caracterise en ce que l'unite (2) 

de determination de la dimension significative de 
l'objet (100) comprend : 

un premier ensemble infrarouge (4, 6) 
constitue d'un premier emetteur infrarouge (4) et d' un 
25 premier recepteur infrarouge (6) , 

un deuxieme ensemble (8, 10) constitue 
d ; un deuxieme emetteur infrarouge (8) et d' un deuxieme 
recepteur infrarouge (10) , 

les deux ensembles infrarouges (4, 6 ; 8, 10) etant 
30 eloignes 1 ' un de 1' autre d'une distance connue (d) , et 
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emettant des faisceaux infrarouges qui sont parralleles 
entre eux, 

et la dimension significative (x) de 
l'objet (100) etant deduite de la reponse inf rarouge 
5 obtenue lorsgue l'objet (100) est deplace de maniere a 
intercepter successivement le premier faisceau. et le 
deuxieme faisceau infrarouges, selon une direction 
sensiblement perpendiculaire a celle des axes (12, 14) 
des deux faisceaux, ladite reponse inf rarouge 
10 correspondantr. a la fraction (24) du deuxieme faisceau 
non encore interceptee par l'objet (100) lorsqu' il 
intercepte encore la moitie (22) du premier faisceau. 

5. Syst^me selon la revendicat ion 4, 
caracterise en ce que 1' unite (2) de determination 

15 comprend en outre un troisieme ensemble (1-6, 18) 
emetteur-recepteur photoelectrique , dispose en a-mont du 
premier ensemble inf rarouge (4, 6) par rapport au 
deuxieme ensemble inf rarouge (8, 10), et destine a 
l'ajustement prealable de l'intensite des deux 

2 0 faisceaux infrarouges. 

6. Systeme selon la revendicat ion 4 ou 5, 
caracterise &zi ce que la dimension significative (x) de 
l'objet (100) est obtenue apres avoir fait cieplacer 
ledit objet QN fois et me sure Q reponses inf rrarouges 

25 RI (q) , pour eg compris entre 1 et Q, par une relation du 
type : 

x = A 4 . (moyenne RI(q)) 4 + A 3 . (moyenne RI(q)) 3 -f 

A 2 . (moyenne RI(q)) 2 + A x . (moyenne RI (q) ) 1 + A 0 , 

ou : 

3 0 P± 0 , Ai," A 2 , A 3 , A 4 sont des coefficients 

obtenus- precedemment en appliquant la meme relation a 
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au moins qiaatre objets etalons en dimension (edim) , 
pour lesquels une reponse infrarouge RI (edim) est 
mesuree . 

7- Systeme selon 1'une quelconque des 
5 revendi cat ions 1 a 6, caracterise en ce que l'urxite 
(30) de determination de 1'intensite attenuee dL'un 
faisceau photohique est une unite de deter mihat ion ' par" 
spectrometry e gamma, qui comprend : 

un ensemble (32) forme d'une source et 
10 d'un collima.teur , 

un ensemble (40) forme d'un detecteuir et 
d ! un collimateur, 

une chaine d' acquisition et de comptage 
des photons gamma (4 8) . 
15 8. Systeme selon la revendication 7, 

caracterise en ce que la chaine d' acquisition^ et de 
comptage (43) comprend : 

un detecteur a germanium de haute 

densite, 

2 0 - un preamplif icateur (50) , 

un module DSP (Digital Signal Processor) 

(52) , 

un module haute tension (54) , 
un module reseau (56) , 
25 - un ordinateur d' acquisition de donnees 

(170), 

un cryostat (60) • 
9. Systeme selon l'une quelconque des 
revendications 1 a 8, caracterise en ce que les moyens 

3 0 de transport (70, 72, 80, 82, 84, 86, 88) comprennent 
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un plateau tournant (70) et un moteur pas a pas (72) 
d'entrainement dudit plateau (70) . 

1 0 . S ys t erne s e 1 on 1 ' une que 1 conque de s 
revendlcations 1 a 9, caracterise en ce que les moyens 

5 de transport comprennent u_n bras de manutention (80) . 

11. SystSme selon la revendication 10, 
caracterise en ce que le bras de manutention (80) est 
un bras articule equipe d/une pince d'extremite (82) 
destinee a saisir et depos er l'objet (100). 

10 12. Syst^me selon l'une quelconque des 

revendlcations 3 a 11, caracterise en ce que les 
premiers moyens d' ajustemexit comprennent : 

- une glissidr-e (74) pour fixer la position 
d'un socle (26) de l'unite (2) de determination de la 

15 dimension significative de l'objet selon une direction 
X. 

- un verin (7 6) pour deplacer le premier 
ensemble infrarouge (4, 6) par rapport le deux i erne 
ensemble infrarouge (8, 10) de ladite unite (2) selon 

2 0 une direction Y perpendiciilaire a la direction Y. 

- un verin (78) pour deplacer ledit socle 
(26) de ladite unite (2) selon une direction Z 
perpendiculaire au plan (XI, Y) . 

13. Systeme selon l'une quelconque des 
2 5 revendlcations 1 a 12, caracterise en ce que les 
deuxi ernes moyens d' ajustement comprennent un support 
d' irradiation (90) sur ILequel est installe l'objet 
(100) entre une source (32) et un detecteur (40) de 
l'unite (30) de determination de l'intensite attenuee 
30 du faisceau photonique traversant l'objet (100). 
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14. Systeme selon la revendication 13, 
caracterise en ce que les deuxi<Smes moyens d'ajustement 
comprennent : 

- une glissiere (94> pour fixer un support 
5 d' irradiation (90) selon une direction X, 

- un verin (96) poxir deplacer, selon une 
direction Y perpendiculaire §. la direction X, ledit 
support d' irradiation (90) entre une source (32) et un 
detecteur (40) de 1' unite (3D) de determination de 

10 l'intensite attenuee du faiscea_Ti photonique traversant 
l'objet (100) , 

- un verin (98) pomr deplacer, selon une 
direction Z perpendiculaire a~u plan (X, Y) , ledit 
support d' irradiation (90) entre une source (32) et un 

15 detecteur (40) de 1' unite (30) de determination de 
l'intensite attenuee du faiscea_Ti photonique traversant 
l'objet (100) . 

15 . Procede de mis e en ceuvre du systeme 
pour la determination automat iqi_ae de la masse volumique 

2 0 d' un objet (10 0) appartenant a un lot d'objets selon 

1'une quelconque des revendic at ions 1 a 14, ledit 
systeme comportant une unite (2) de determination d'une 
dimension significative (x) d' un objet (100) et une 
unite (300) de determination c3.e l'intensite (I) d'un 
25 faisceau photonique attenue par la traversee dudit 
objet (100), caracterise en ce qu'il comprend les 
etapes d'etalonnage suivantes : 

- une etape 1 d'etalonnage de position de 
deux ensembles inf rarouges- (4, 6 ; 8, 10) de 1 ' unite 

3 0 (2) de determination de la dimension significative des 

objets (100) , 



WO 2005/100952 



PCT/FR2005/000838 



53 



- une etape 2 d' etalonnage de position d'un 
support d' irradiation (90) de 1.' unite (3 0) de 
determination par spectrom^trie gamma de 1'intensite du 
faisceau photonique attenuee par la traversee des 

5 objets (100) , 

- une etape 3 d'etalonnage de mesure de 
1' ensemble source -detecteur (32, 40) c3.e 1 ' unite (30), 

et en ce qu'il compread des etapes de 
determination proprement dite d& la dimension 
10 significative (x) des objets (10 0) , qxai sont effectuees 
sur chaque objet (100) dudit lot d' objets. 

16. Procede selon la revendication 15, 
caracterise en ce que 1 ' etape 1 d' etalonnage comporte 
une saisie par l'operateur d'un ensemble de parametres 
15 d' entree au moyen d'un module interactif, parmi 
lesquels : 

- la configuration de composants qui ont un 
deplacement micrometrique, parmi lesquels deux verins 
(76, 78), 

2 0 la configuration du plateau tournant 

(70), c'est-a-dire la nature des objets qui occupent 
differents emplacements prevus sur celui-ci, 

1 1 emplacement occupe par chaque objet 
etalon en dimension (edim) sur le plateau tournant 
25 (70), 

- une position (Z meS ure) d'un socle (26) de 
1' unite (2) suivant la direction Z, 

- des positions (Y(l), Y (N) ) bornant un 
intervalle de deplacement du premier ensemble 

3 0 infrarouge (4, 6) suivant la direction Y, 
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- un pas (INT) de displacement du premier 
ensemble infrarouge (4, 6) suivant la direction Y. 

17. Procede selon la revendication 16, 
caracterise en ce que l'etape 1 d'etalonnage comporte 

5 en outre les operations suivantes : 

a) deplacement du socle (26) suivant la 
direction Z jusqu'a la position (Z meS ure) 

b) deplacement angulaire du plateau 
tournant (70) afin de transporter un objet etalon en 

10 dimension (edim) jusqu'a sa position initiale mesure 
par rapport a l'unite (2), 

c) deplacement du premier ensemble 
infrarouge (4, 6) suivant la direction Y jusqu'a sa 
position de depart (Y(l)), 

15 d) deplacement progressif, par increments 

successifs de (INT) , du premier ensemble infrarouge (4, 
6) suivant la direction Y 7 en l'eloignant du deuxieme 
ensemble infrarouge (8, 10) fixe a une position Y FIX/ 
entre les positions Y(l) et Y (N) , et determination 

20 simultanee de la reponse infrarouge (RI (n) ) de 1' objet 
(edim), correspondant a chaque position (Y (n) ) 

e) calcul de la reponse infrarouge optimale 

RIopt, 

f) calcul de la position optimale Yqpt du 
25 premier ensemble infrarouge (4, 6) par rapport au 

deuxieme ensemble infrarouge (8, 10). 

18. Procede selon la revendication 17, 
caracterise en. ce que 1' operation d) de deplacement 
progressif comporte les sous -operations suivantes : 

3 0 d-1) deplacement angulaizre du 

plateau tournant (70) afin de transporter 1' objet 



WO 2005/100952 



PCT/FR2005/000838 



55 



etalon en dimension (edim) a sa position f inale de 
me sure, 

d-2) mesure de la reponse infrarouge 
(RI (n) ) dudit objet etalon en dimension (edim), 
5 d-3) deplacement angula-ire du 

plateau tournant (70) afin de ramener 1' objet etalon en 
dimension (edim) a sa position initiale de mesixre. 

19. Procede selon la revendication 17 ou 
18, caracterise en ce que la reponse infrarouge 

10 optimale RIopt est obtenue par la relation : 

^opt 2 

ou : RImin est la valeur de la saturation minimale de la 

reponse infrarouge , 
et : RImax est la valeur de la saturation maximal e de la 
15 reponse infrarouge. 

20. Procede selon l'une quelcoinque des 
revendi cat ions 17 a 19, caracterise en ce que 
1 ' operation f) de calcul de la position opti_male Y 0 pt 
est obtenue de la maniere suivante : 

RI OPT -RIG) 

2 0 si <1 , alors Y 0 pt = Y(j) 

RI OFr -RI(k) 

sx >1 , alors Yopt = Y(k) 

RI OFr -RI(k) 

ou RI ( j ) et RI (k) sont deux valeurs prece demment 
calculees de la reponse infrarouge entre lesguelles est 
comprise la reponse optimale recherchee RZEqpt, qui 
25 correspondent respect ivement a deux positions Y(j) et 
Y(k) du premier ensemble infrarouge (4, 6). 

21. Procede selon l'une quelcoinque des 
revendi cat ions 15 a 20, caracterise en ce que H'etape 2 
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d'etalonnage comporte une saisie par 1 7 operateur d'un 
ensemble de param^tres d' entree au moyen d'un module 
interactif, parmi lesquels : 

- la configuration de composants qui ont un 
5 deplacement micrometrique, parmi lesquels deux verins 

(96, 98), 

la configuration du plateau tournant 
(70), c ' est-a-dire la nature des objets qui occupent 
differents emplacements sur celui-ci, 
10 - 1' emplacement occupe par chaque objet 

etalon en masse volumique (emas) sur le plateau 
tournant (70) , 

- une duree de mesure ou temps de comptage, 

des positions Z(l) et Z (N) bornant un 
15 intervalle de deplacement d'un support d' irradiation 
(90) suivant la direction Z, 

- un nombre M de mesures de l'intensite 
photonique attenuee par la traversee de chaque objet 
etalon en masse volumique, pour chaque position Z(i) 

20 occupee par le support d' irradiation, pour i = 1,...,N. 

22. Procede selon la revendication 21, 
caracterise en ce que 1 7 etape 2 d'etalonnage comporte 
en outre les operations suivantes : 

a) determination de la dimension 
25 significative (x mas ) de chaque objet etalon en masse 

volumique (emas) , 

b) deplacement angulaire du plateau 
tournant (70) , d'un angle (A) , afin de transporter 
ledit objet etalon en masse volumique (emas) en une 

3 0 position intermediaire ou il est saisi par un bras de 
prehension (8 0) , 
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c) positionnement dudit objet (emas) sur un 
support d' irradiation (90), 

d) ajustement proprement dit de la position 
du support d' irradiation (90) par rapport a une source 

5 (32) et un detecteur associe (40) , 

e) transport retour de 1' objet etalon en 
masse volumique (emas) sur le plateau tournant (70) , 
par une operation inverse a 1 7 operation c) . 

23. Procede selon la revendications 22, 
10 caracterise en ce que 1' operation c) de positionnement 
de 1 ' objet (emas) sur le support d' irradiation (90) 
comporte les sous-operations suivantes : 

c-1) deplacement du support d' irradiation 
(90) , suivant la direction Z et vers le bas, 
15 c-2) deplacement du bras de prehension (80) 

depuis sa position d' attente jusqu'a la verticale de la 
position intermediaire de l'objet (emas), 

c-3) saisie, par le bras de prehension 
(80), de l'objet emas, et transport de celui-ci jusqu'a 
2 0 la verticale d'une face superieure (92) du support 
d' irradiation (90) , 

c-4) deplacement du support d' irradiation 
90 jusqu'a la position Z(l), suivant la direction Z et 
vers le haut , 

25 c~5) depose de l'objet (emas) sur la face 

superieure (92) du support d' irradiation (90), par le 

bras de manutention 80, 

c-6) deplacement retour du bras de 

manutention (80) jusqu'a sa position d'attente, 
30 c-7) mise en butee, suivant la direction Y, 

de l'objet (emas) sur la face superieure (92). 
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24. Procede selon la revendi cat ions 22 ou 
23 , caracterise en ce que 1 ' operation d) d' ajustement 
proprement dit de la position du support d' irradiation 
(90) par rapport a la source (32) et au detecteur 

5 associe (40) comporte les sous -operations suivantes : 

d-1) deplacement progressif du support 
d' irradiation (90) suivant la direction Z entre deux 
positions Z(l) et Z (N) predeterminees , 

d-2) pour chaque position Z(i), i = 1, 
10 N, irradiation de l'objet etalon en masse volumique 
(emas) par le faisceau photonique, un nombre M de fois, 
ce qui conduit a un ensemble de valeurs d'intensite 
attenuee I(i, j), ou i = 1, N represente le nombre 

de positions successives Z(i) occupees par le support 
15 d' irradiation (90) et j = 1, M represente le nombre 

d' irradiations effectuees a chaque position Z(i), 

d-3) calcul de la position optimale (Z 0 pt) 
du support d' irradiation (90) a partir d'une regression 
polynomiale d'ordre 4 des positions Z(i) par rapport 
20 aux intensites attenuees I(i, j), cette regression 
polynomiale d'ordre 4 etant predeterminee et integree 
comme une donnee d'une unite (200) d' acquisition, de 
traitement et d' analyse. 

25. Procede selon 1'une quelconque des 
25 revendications 15 a 24, caracterise en ce que l'etape 3 

d'etalonnage de 1'unite 30 de determination par 
spectrometrie gamme comprend les operations 
automatisees suivantes : 

a) mesure de l'intensite photonique (Iemas) 
3 0. attenuee par la traversee d'un objet etalon en masse 
volumique (emas) 
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b) calcul du coefficient massique 
ci ' attenuation (|Ll m ) de l'objet etalon en masse 

volumique, par la relation : = L 1 



H'm X emas ^0 

26. Procede selon l'une quelconque des 
5 revendications 15 a 25, caracterise en ce que les 

etapes de determination proprement dite comprennent : 

une etape 4 de determination de la 
dimension significative (x) de l'objet (100) a 
contrSler , 

10 - une etape 5 de transport de l'objet 

(100) vers un support d' irradiation (90), 

une etape 6 d'ajustement de la position 
de l'objet (10 0) par ajustement de la position du 
support d' irradiation (90) par rapport a une source 
15 (32) et un detecteur associe (40) , 

une etape 7 de determination de 
l'intensite attenuee (I) du faisceau photonique 
transmis a travers l'objet (100) , 

une etape 8 d' acquisition, traitement et 
2 0 analyse du spectre obtenu, 

une etape 9 de determination de la 

variation relative de masse volumique (p) de 

P 

l'objet (100) par rapport a celle d'un ou plusieurs 
objet(s) etalon(s) en masse volumique (emas) , 
25 - une etape 10 de transport ret our de 

l'objet (100) jusqu'a son emplacement sur le plateau 
tournant (70) . 

27. Procede selon la revendication 26, 
caracterise en ce que 1' etape 4 de determination de la 
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controler comporte une saisie par l'operateur d'un 

ensemble de parametres d' entree au moyen d'un module 

interactif, parmi lesquels : 
5 la configuration d'un plateau tournant 

(70), c'est-a-dire la nature des objets qui occupent 

differents emplacements sur celui-ci, 

1 ' emplacement occupe par l'objet (100) 

sur le plateau tournant (70) , 
10 - une position (Z meS ure) suivant la direction 

Z du socle (26) de 1'unite (2), 

- un nombre P de mesures infrarouges pour 

chaque objet etalon en dimension (edim(n)), n = 1, 

N, ou N est le nombre d' objets etalons en dimension, 
15 - un nombre Q de mesures infrarouges pour 

l'objet (100) . 

28. Precede selon la revendication 27, 

caracterise en ce que l'etape 4 de determination de la 

dimension significative (x) de l'objet (100) a 
2 0 controler comporte en outre" les operations suivantes : 

a) deplacement d ; un socle (26) de 1' unite 
(2) suivant la direction Z jusqu'a la position (Z mesure ) , 

b) deplacement d'un premier ensemble 
infrarouge (4, 6) suivant la direction Y, jusqu'a une 

25 position (Y mesU re) definie par : Y mesU re = Y 0 pt + (x e dim - 

XeditnMOY) / OU : 

Yopt est la position optimale provenant de 
l'etape 1 d' etalonnage , 

Xedim est la dimension de l'objet etalon en 
■30 dimension (edim) utilise lors de l'etape 1 
d' etalonnage , 
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XedimMOY est la dimension significative 
moyenne de tous les objets (edim) gtalons en dimension, 

c) mesure repetee P fois de la reponse 
infrarouge RI (p) , p = 1, P des N objets etalpns en 

5 dimension (edim(n)), n = 1, ... N, ce qui conduit a un 
ensemble de valeurs RI (n, p) , 

d) calcul proprement dit de la dimension 
significative de l'objet (100). 

29. Procede selon la revendication 28, 
10 caracterise en ce que 1' operation d) de calcul 

proprement dit de la dimension significative (x) de 
l'objet (100) est effectuee de la maniere suivante : 

d-1) raise en oeuvre d'une regression 
polynomiale d'ordre 4 des dimensions signif icatives 
15 Xedim(n) de chacun des N objets etalons en dimension 
edim(n) , en fonction de la moyenne des reponses 

V RI(n,p) 

infrarouges RI edimMOY = — — dudit objet etalon en 

dimension edim(n) , pour calculer les coefficients A 0/ 
Ai, A 2 , A 3/ A 4 d'une relation du type : 

20 Xedim(n) = A 4 . (RIedimMOY (n) ) 4 + A 3 . (RIedimMOY (n) ) 3 + 

A 2 . (RIedimMOY (n) ) 2 + A x . (RIedimMOY (n) ) 1 + A 0 , 

d-2) mesure repetee Q fois de la reponse 
infrarouge RI (q) , q = 1, Q de l'objet a controler 

T Ri(q) 

100, calcul de la moyenne RI- — de ces reponses 

25 infrarouges, et calcul de la dimension significative x 
de l'objet 10 0 par la relation : 

x = A 4 . (RI) 4 + A 3 . (RI) 3 + A 2 . (RI) 2 + Ai . (RI) 1 + A 0 

30. Utilisation du systeme selon l'une 
quelconque des revendications 1 a 14 et du procede 
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selon l'une quelconque des revendi cat ions 15 a 29, pour 
le contrdle d'objets (100) en cours de fabrication. 

31. Utilisation selon la revendication 30, 
dans laquelle les objets (100) sont des pastilles de 
5 combustible nucleaire. 
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FIG. 8 
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ETALONNAGES PREALABLES 
Etape 1 

etalonnage de position des deux ensembles infrarouges de 
I'unite de determination de la dimension significative de I'objet 

Etape 2 

etalonnage de position du support d'irradiation de I'unite 
de determination de Pintensite du faisceau photonique attenuee 
par la traversee de I'objet 

Etape 3 

etalonnage de mesure de I'ensemble source-detecteur de I'unite 
de determination de I'intensite du faisceau photonique attenuee 
par la traversee de I'objet 



FIG. 9A 
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DETERMINATION PROPREMENT DITE 

DE LA VARIATION RELATIVE DE MASSE VOLUMIQUE 

Etape 4 

determination de la dimension significative de I'objet a contrdler 
Etape 5 

transport de I'objet vers le support d'irradiation 
Etape 6 

ajustement de la position de I'objet par ajustement de la position 
du support d'irradiation par rapport a une source et un detecteur 
associe 

Etape 7 

determination de I'intensite attenuee du faisceau 
photonique transmis a travers I'objet 

Etape 8 

acquisition, traitement et analyse du spectre obtenu 
Etape 9 

determination de la variation relative de masse volumique 
de I'objet par rapport a celle d'un ou plusieurs objet(s) 
etalon(s) en masse volumique 

Etape 1 0 

transport retour de I'objet jusqu'a son emplacement sur 
le plateau tournant 



FIG. 9B 
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Etape 1 



saisie manuelle des parametres d'entree 



4 

sequence cooperations automatisees : 



a) deplacement du socle suivant la direction Z 
jusqu'a la position Z mesU re 

b) deplacement angulaire du plateau tournant 
pour amener I'objet edim a sa position initiale 
de mesure 

c) deplacement du premier ensemble infrarouge 
suivant la direction Y jusqu'a sa position de 
depart Y(1) 

d) deplacement du premier ensemble infrarouge 
suivant la direction Y, par increments de 

INT, et determination simultanee de la reponse 
infrarouge Rl(n) de I'objet edim, correspondant 
a chaque position Y(n), pour n = 1, N 



e) calcul de la reponse infrarouge optimale 

pi RImax~R*min 
ki opt 



f) calcul de la position optimale Yqpt du 
premier ensemble infrarouge par rapport au 
deuxieme ensemble infrarouge, et de 
I'ecartement optimal d entre les deux ensembles 
infrarouges 



FIG. 10 
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Etape 2 



saisie manuelle des 



parametres d 'entree 



v 



sequence d'operations automatisees : 

a) mesure de la dimension significative x emas de 
I'objet etalon en masse volumique 

b) deplacement angulaire du plateau tournant, 
pour amener I'objet etalon en masse volumique 
emas en une position intermediate ou ii est 
saisi par le bras de prehension, 

c) positionnement de Pobjet emas sur le support 
d 'irradiation 

d) ajustement proprement dit de la position du 
support d'irradiation par rapport a I'ensemble 
source-detecteur : 

d-1) deplacement progressif du support 
d'irradiation suivant la direction Z entre les positions 
Z(1) et Z(N) predetermines 

d-2) pour chaque position Z(i), 

irradiation M fois de I'objet emas par le 
faisceau photonique, et obtention des valeurs 
d'intensite attenuee l(i, j) 

i = 1, N = nombre de positions Z(i) entre 
Z(1) et Z(N) 

j = 1 M a nombre d'irradiations pour chaque 

position Z(i) 

d-3) calcul de la position optimale Z opt 
du support d'irradiation a partir des positions 
Z(i) et des intensites attenuees l(i, j) 

e) transport retour de I'objet emas sur le 
plateau tournant 




11 
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Etape 3 



a) mesure de I'intensite photonique l emas 
attenuee par la traversee d'un objet etalon en 
masse volumique emas pris comme reference 



b) calcul du coefficient massique d'attenuation u. m 
de I'objet etalon en masse volumique par la relation 

0 \ t ^emas 

remas LI X " n t 

emas J- q 



FIG. 12 
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Etape 4 



saisie manuelles des parametres d'entree 



± 

sequence d'operations automatisees : 

a) deplacement du socle suivant la direction Z 
jusqu'a la position Z mesure 

b) deplacement du premier ensemble infrarouge 
suivant la direction Y, jusqu'a la position 

Xnesure ~ Y OPT + ( x edim ~ x edimMOY) 

c) mesure repetee P fois de la reponse 
infrarouge Rl(p), p = 1, P des N objets 
etalons en dimension edim(n), n = 1, ... N, 
conduisant a un ensemble de valeurs Rl(n, p), 

d) calcul proprement dit de la dimension 
significative x de I'objet 100 : 

d-1) calcul de la moyenne Rl edimMoydes P reponses 
infrarouges de chacun desdits objets edim(n) 
dont les dimensions significatives x edim (n) sont 
connues, puis calcul des coefficients A o , A 1r 
A 2 , A 3, A4 de la relation 

Xedim 0"0 = A 4 . ( Rl edimMOY (n)) 4 + A3 . (Rl edimMOY 
(n)) 3 + A 2 . (Rl edimMOY (n)) 2 + A1 . (Rl edimMOY (n)) 1 
+ A 0 

d-2) mesure de la reponse infrarouge Rl(q), q = 
1, Q de I'objet 100, puis calcul de la 
moyenne des Rl(q), et calcul de la dimension 
significative x de I'objet 100 par la relation : 

x = A 4 . (Rl) 4 + A 3 . (Rl) 3 + A 2 . (Rl) 2 + Ai . 
(Rl) 1 + A 0 . 



FIG. 13 
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Etape 9 



calcul automatise de la variation relative 

de masse volumique de I'objet 100 par rapport a 
celle d'un ou plusieurs objet(s) etalon(s) en 
masse volumique emas a partir de la relation: 



P : 
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